Multi-Scale Modeling of Hypertension

The focus of this work is the coupling of a 1-D lumped parameter model representing the circulatory system to a 3-D finite element based left ventricle (LV) model in order to study the effects of mild hypertension on the cardiovascular system. A Finite Element LV model under normotensive loading (116/80 mmHg) was developed as well as a mild hypertension (165/90 mmHg) model. In both cases, coupled analysis was utilized so that at one diastolic time point and four systolic time points the values for the LV volumes and pressures were determined from the steady-state JSim solution. The normotensive model had an average first principal stress of 39.1KPa while the hypertensive case showed an increased value of 51.8KPa representing a 32.3 % increase. A relatively mild increase in the afterload resulted in a pronounced increase in workload to maintain the same systemic flow. 1

Graphical Approach to Model Reduction for Nonlinear Biochemical Networks

Model reduction is a central challenge to the development and analysis of multiscale physiology models. Advances in model reduction are needed not only for computational feasibility but also for obtaining conceptual insights from complex systems. Here, we introduce an intuitive graphical approach to model reduction based on phase plane analysis. Timescale separation is identified by the degree of hysteresis observed in phase-loops, which guides a “concentration-clamp” procedure for estimating explicit algebraic relationships between species equilibrating on fast timescales. The primary advantages of this approach over Jacobian-based timescale decomposition are that: 1) it incorporates nonlinear system dynamics, and 2) it can be easily visualized, even directly from experimental data. We tested this graphical model reduction approach using a 25-variable model of cardiac β1-adrenergic signaling, obtaining 6- and 4-variable reduced models that retain good predictive capabilities even in response to new perturbations. These 6 signaling species appear to be optimal “kinetic biomarkers” of the overall β1-adrenergic pathway. The 6-variable reduced model is well suited for integration into multiscale models of heart function, and more generally, this graphical model reduction approach is readily applicable to a variety of other complex biological systems

Multifractality in Human Heartbeat Dynamics

Recent evidence suggests that physiological signals under healthy conditions may have a fractal temporal structure. We investigate the possibility that time series generated by certain physiological control systems may be members of a special class of complex processes, termed multifractal, which require a large number of exponents to characterize their scaling properties. We report on evidence for multifractality in a biological dynamical system --- the healthy human heartbeat. Further, we show that the multifractal character and nonlinear properties of the healthy heart rate are encoded in the Fourier phases. We uncover a loss of multifractality for a life-threatening condition, congestive heart failure.Comment: 19 pages, latex2e using rotate and epsf, with 5 ps figures; to appear in Nature, 3 June, 199

Nonlinear complexity analysis of brain fMRI signals in schizophrenia

Peer reviewedPublisher PD

CytoSolve: A Scalable Computational Method for Dynamic Integration of Multiple Molecular Pathway Models

A grand challenge of computational systems biology is to create a molecular pathway model of the whole cell. Current approaches involve merging smaller molecular pathway models’ source codes to create a large monolithic model (computer program) that runs on a single computer. Such a larger model is difficult, if not impossible, to maintain given ongoing updates to the source codes of the smaller models. This paper describes a new system called CytoSolve that dynamically integrates computations of smaller models that can run in parallel across different machines without the need to merge the source codes of the individual models. This approach is demonstrated on the classic Epidermal Growth Factor Receptor (EGFR) model of Kholodenko. The EGFR model is split into four smaller models and each smaller model is distributed on a different machine. Results from four smaller models are dynamically integrated to generate identical results to the monolithic EGFR model running on a single machine. The overhead for parallel and dynamic computation is approximately twice that of a monolithic model running on a single machine. The CytoSolve approach provides a scalable method since smaller models may reside on any computer worldwide, where the source code of each model can be independently maintained and updated

Bolus characteristics based on Magnetic Resonance Angiography

BACKGROUND: A detailed contrast bolus propagation model is essential for optimizing bolus-chasing Computed Tomography Angiography (CTA). Bolus characteristics were studied using bolus-timing datasets from Magnetic Resonance Angiography (MRA) for adaptive controller design and validation. METHODS: MRA bolus-timing datasets of the aorta in thirty patients were analyzed by a program developed with MATLAB. Bolus characteristics, such as peak position, dispersion and bolus velocity, were studied. The bolus profile was fit to a convolution function, which would serve as a mathematical model of bolus propagation in future controller design. RESULTS: The maximum speed of the bolus in the aorta ranged from 5–13 cm/s and the dwell time ranged from 7–13 seconds. Bolus characteristics were well described by the proposed propagation model, which included the exact functional relationships between the parameters and aortic location. CONCLUSION: The convolution function describes bolus dynamics reasonably well and could be used to implement the adaptive controller design

Relational grounding facilitates development of scientifically useful multiscale models

We review grounding issues that influence the scientific usefulness of any biomedical multiscale model (MSM). Groundings are the collection of units, dimensions, and/or objects to which a variable or model constituent refers. To date, models that primarily use continuous mathematics rely heavily on absolute grounding, whereas those that primarily use discrete software paradigms (e.g., object-oriented, agent-based, actor) typically employ relational grounding. We review grounding issues and identify strategies to address them. We maintain that grounding issues should be addressed at the start of any MSM project and should be reevaluated throughout the model development process. We make the following points. Grounding decisions influence model flexibility, adaptability, and thus reusability. Grounding choices should be influenced by measures, uncertainty, system information, and the nature of available validation data. Absolute grounding complicates the process of combining models to form larger models unless all are grounded absolutely. Relational grounding facilitates referent knowledge embodiment within computational mechanisms but requires separate model-to-referent mappings. Absolute grounding can simplify integration by forcing common units and, hence, a common integration target, but context change may require model reengineering. Relational grounding enables synthesis of large, composite (multi-module) models that can be robust to context changes. Because biological components have varying degrees of autonomy, corresponding components in MSMs need to do the same. Relational grounding facilitates achieving such autonomy. Biomimetic analogues designed to facilitate translational research and development must have long lifecycles. Exploring mechanisms of normal-to-disease transition requires model components that are grounded relationally. Multi-paradigm modeling requires both hyperspatial and relational grounding

Comparison of four mathematical models to analyze indicator-dilution curves in the coronary circulation

While several models have proven to result in accurate estimations when measuring cardiac output using indicator dilution, the mono-exponential model has primarily been chosen for deriving coronary blood/plasma volume. In this study, we compared four models to derive coronary plasma volume using indicator dilution; the mono-exponential, power-law, gamma-variate, and local density random walk (LDRW) model. In anesthetized goats (N = 14), we determined the distribution volume of high molecular weight (2,000 kDa) dextrans. A bolus injection (1.0 ml, 0.65 mg/ml) was given intracoronary and coronary venous blood samples were taken every 0.5–1.0 s; outflow curves were analyzed using the four aforementioned models. Measurements were done at baseline and during adenosine infusion. Absolute coronary plasma volume estimates varied by ~25% between models, while the relative volume increase during adenosine infusion was similar for all models. The gamma-variate, LDRW, and mono-exponential model resulted in volumes corresponding with literature, whereas the power-model seemed to overestimate the coronary plasma volume. The gamma-variate and LDRW model appear to be suitable alternative models to the mono-exponential model to analyze coronary indicator-dilution curves, particularly since these models are minimally influenced by outliers and do not depend on data of the descending slope of the curve only

In silico assessment of biomedical products: the conundrum of rare but not so rare events in two case studies

In silico clinical trials, defined as “The use of individualized computer simulation in the development or regulatory evaluation of a medicinal product, medical device, or medical intervention,” have been proposed as a possible strategy to reduce the regulatory costs of innovation and the time to market for biomedical products. We review some of the the literature on this topic, focusing in particular on those applications where the current practice is recognized as inadequate, as for example, the detection of unexpected severe adverse events too rare to be detected in a clinical trial, but still likely enough to be of concern. We then describe with more details two case studies, two successful applications of in silico clinical trial approaches, one relative to the University of Virginia/Padova simulator that the Food and Drug Administration has accepted as possible replacement for animal testing in the preclinical assessment of artificial pancreas technologies, and the second, an investigation of the probability of cardiac lead fracture, where a Bayesian network was used to combine in vivo and in silico observations, suggesting a whole new strategy of in silico-augmented clinical trials, to be used to increase the numerosity where recruitment is impossible, or to explore patients’ phenotypes that are unlikely to appear in the trial cohort, but are still frequent enough to be of concern